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ABSTRACT 

Context. Triggered star formation has been discussed for many years, and evidence for the formation of stars and cores 
triggered by Hn regions is under debate. 

Aims. We investigate the imprints of an expanding Hn region on a pre-existing starless clumps. 

Methods. We selected an infrared dark filament spanning 0.8°. One portion of this filament, G18.93-0.03 is a prominent 
dust complex, with the molecular clump G18.93/m being IR dark from near IR wavelength up to 160 /im. Spitzer 
composite images show an IR bubble spatially associated with G18.93-0.03. We use GRS 13 CO and IRAM 30m H 13 CO + 
data to disentangle the large and small scale spatial structure of the region. From ATLASGAL submm data we calculate 
the gas mass, while we use the H 13 CO + line width to estimate its virial mass. To study the IR properties of G18.93/m we 
use HERSCHEL data and produce temperature maps from fitting the spectral energy distribution. With the MAGPIS 
20 cm and SuperCOSMOS H Q data we trace the ionized gas, and the VGPS Hi survey provides information on the 
atomic hydrogen gas. 

Results. We show that the bubble is spatially associated with G18.93, located at a kinematic near distance of 3.6 kpc. 
The total gas mass of ~ 870 M splits up into 6 sub-clumps, of which G18.93/m is the most massive with 280 M©. The 
virial analysis shows that it may be gravitationally bound and has neither Spitzer young stellar objects nor mid-IR 
point sources within. Therefore we call it pre-stellar. Fitting the spectral energy distribution reveals a temperature 
distribution that decreases towards its center, but heating from the ionizing source puts it above the general ISM 
temperature. We find that the bubble is filled by Hn gas, ionized by an 08.5 star. Between the ionizing source and 
the IR dark clump G18.93/m we find a layered structure of hydrogen phases, from ionized to atomic to molecular gas, 
revealing a photon dominated region. Furthermore, we identify an additional velocity component within the bubble's 
8 /im emission rim at the edge of the infrared dark cloud and speculate that it might be shock induced by the expanding 
Hll region. 

Conclusions. While the elevated temperature allows for the build-up of larger fragments, and the shock induced velocity 
component may lead to additional turbulent support, the density profile of G18.93/m does not show signatures of the 
expanding bubble. While the first two signatures favor high-mass star formation, we do not find conclusive evidence 
that the massive clump G18.93/m is prone to collapse because of the expanding Hn region. 
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1. Introduction 



Based on th e Galactic plane survey GLIMPSE ( |Benjamiri| 



2003), it has become possible to study bubble-like 



|et al 

structures in the Milky W ay in a statistical sense. Initially, 
(2006k found a 'Bubbling Galactic Disk' 
ubbles. Today, several compilations of 



Churchwell et al. 
and 



identified 322 

bubbles in the Milky Way exist (e.g Churchwell et al. , 2006| 
2007). The most extensive of which has been co mpiled by 



2012), a citizen 



7 THe Milky Way Project' ( [Simpson et al. 
science project. They visually identified ancl classified more 
than 5000 bubbles. Infrared bubbles are defined by their 
bright 8 /im emission rims. Although their origin may be 
manifold, many of the bubbles have shown to be created by 
a, not necessarily central, ionizing source. In this context, 



the emission rim originates from UV excited polycyclic aro- 
matic hydrocarbons (PAHs). While all hydrogen ionizing 
photons are absorbed within the borders of an ionization 
front, lower energy UV photons pass through the ionization 
front and may excite the PAHs. In addition, the ionizing 
source heats its surrounding gas, which expands and drives 
a shock front beyond the ionization front. In between the 
shock and ionization front, neutral material can be accu- 



mulated. Elmegreen fc Lada|(|1977|) discussed the triggering 



of star formation on the borders of expanding Hll regions. 
Among others, 'collect and collapse' describes the collapse 
of swept up mat erial along a compress ed layer. |Deharv eng 



et al. 



2003) and Zavagno et al. (2006) established 'collect 
and collapse' observation ally. A~brief overview on triggered 
star formation is given in Deharveng et al. (2010). 
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Other results of the GLIMPSE survey are systematical 
and sensitive studies of infrared dark clouds (IRDCs, e.g. 
|Peretto fc Fuller 2009). Previously discovered by the two 
infrared satellites ISO and MSX, their high column densi- 
ties absorb the background emission even at infrared wave- 
lengths and thus become visible as dark patches. With such 
high column densities, IRDCs are believed to be the cra- 
dles of the next gen eration of stars (|Rathborne et al.] |2006| 



2.1. The cold dust tracing the molecular hydrogen, 
ATLASGAL 

Molecular hydrogen, H2, is hard to trace directly, and no 
pure rotational transitions are observable from the ground. 
However, in regions of high particle densities the gas ther- 
mally couples with the dust. Cold dust can be observed 
by its thermal radiation in the Rayleigh- Jeans regime at 
submm and mm wavelength. The APEX Telescope Large 



Simon et al. 2006). In addition, a rising number ofTarge 
scale Galactic Plane surveys allow one to trace the various 

stages of star formation. From mm and sub-mm observa- |2009[) covers the full inner Galactic plane at 870 fim with 



2009 



tions of cold , dense clumps and cores ( Schuller et al. 
Bally et al. , 2010[) to mid- to near-IR observations of youn^ 



Area Survey of the Galaxy (ATLASGAL, Schuller et al. 



stellaF^FJects 



eniamin et al. , 2003 



and stars as we 



Carey et al. , 2009 



1 as the ionized anc 



Molinari et al. t 

molecular gas (Condon et al.[ |1998 Jackson et al. 2006 



Purcell fc Hoare[ |2010p , an almost complete picture can be 
obtained for regions within the Galactic Plane. 

However, although both filamentary structures and 
IRDCs have gotten much interest over the last years, ex- 



treme IRDCs like the 80 pc long "Nessie" (Jackson et al. 



2010) are still rare objects. A similar large scale filament o: 



more t han 5 pc length has been presented by Kainulainen 



et al. (2011). Figure [T] shows this filament, spanning over 



> 0.87°. At a kinematic near distance of ~ 3.6_ ; 5 kpc, this 
translates to a leng th of ~ 54 pc. (For a discussion of the dis- 
tance, see Sec. 3.2 ) Its coherent velocity structure is shown 
in the bottom panel of Fig. [I] Following the low density gas 
traced by 13 CO along the extinction patch, the smooth ve- 
locity tran sitions suggest that the e ntire region is spatially 
connected. |Kainulainen et al. (2011) studied the extinction 
map of the region and estimate the total gas reservoir to 
be 4.7 x 1O 4 M . Along the filament at lower longitudes, 
IRDC 18223 is a well-studied prot otype IRDC with on go- 
ing high-mass star formation (e.g. Beuther et al.polo and 
references therein). 

A particularly interesting region in the context of trig- 
gered star formation is the prominent IRDC G18.93. Being 
part of the large scale filament described above, its location 
is also coincident with the rim of a IR bubble just above 
the filament (indicated by the upper yellow ellipse in Fig. 
pi. We will show that its dense gas, the filament and the 
oubble have coherent velocities and therefore are all at the 
same distance. Using CLUMPFIND on 870 /im dust con- 
tinuum data we identify six sub-clumps. The most massive 
clump, denoted G18.93/m, is even infrared dark from near 
infrared wavelengths up to 160 /im. 

For the particular case of G18.93/m we will address the 
following questions: Do we see evidence that the bubble 
influences G18.93/m? Do the properties of G18.93/m differ 
from other high-mass starless clumps? 

We will first introduce the data used for the analysis and 
their implications. In Sec. [3] we introduce G18.93 and the 
neighboring bubble, and supplement a more detailed de- 
scription in Sec. |3.3| and onwards. Sec. [4] describes the signs 
of triggered star formation we find and discusses their con- 
sequences. We summarize the paper in a concluding section. 



2. Description and implications of observations and 
data 

We used various publicly available surveys, as well as IRAM 
30 m follow-up observations of the specific region. 



a resolution of 19.2 // and an rms below 50m Jy. Using the 
Larg e APEX Bolometer CAmera (LABOCA, |Siringo et al.| 
2009), each position of the inner Galactic plane has been 



mapped on-the-fly twice, each in different scanning direc- 
tions. The pointing is on the order of 4"and the calibration 
uncertainty is lower than 15% . A detailed descript ion of 
the data reduction is given in Schuller et al.| ( |2009[ ). The 
submm emission is mostly optically thin. 



2.2. IRAM 30 m observations 

We used the IRAM 30 m telescope to map the region 
around G 18. 9 3 i dentified in the ATLASGAL 87 fim images 



Schuller et al. 2009 Tackenberg et al. 2012) in H 13 CO+ 



and SiO. 

H 13 CO + is a well established and common dense gas 
tracer with a critical density at 20 K of ^1.8x 10 5 cm -3 . 
It is generally optically t hin and yields informati on on the 
bulk motion of the gas (Vasyunina et al. , 2011| ). SiO be- 
comes released from the dust grains due to shocks and is 
therefore a we ll known tracer of sh ocked gas and molecular 
outflows (e.g. [Schilke et aT||1997). 



During our observing run in August 2010 the weather 
conditions were poor with average precipitable water 
vapour (pwv) > 10 mm. Although the weather affected the 
pointing, it is, nonetheless, better than a third of the 3 
mm beam (or ~ 9 // ). We used on-the-fly mapping on 5 
boxes to cover all 870 /im continuum emission above 0.3 Jy, 
or 2 x 10 22 cm -2 (assuming a genera l dust temperature of 
15 K, for further details see Sec. |3. 3.1] see Fig.[2]for the cov- 
erage of the observations). The IRAM 30 m beam width for 
both H 13 CO + and SiO is 29.9 // and the 3 mm setup chosen 
provides a spectral resolution of 40 kHz at 40 MHz band- 
width. This translates to a native velocity resolution of 
0.14km/s at 86 GHz. The data reduction was done using 
CLASS from the GILDASQ package. We subtracted a lin- 
ear baseline and removed the three central channels to avoid 
occasional spikes from the backends. The maps have been 
produced with a pixel scale of 14.2". With a total integra- 
tion time of 11.3 h the average rms at a velocity resolution 
of 0.25 km/s is 0.08 K and 0.05 K for H 13 CO+ and SiO, re- 
spectively. 

Simultaneously with the 3 mm observations, we 
recorded two of the H 2 CO(3-2) lines, CH 3 OH(4-3), and 
HCsN(24-23) at 1 mm. However, due to the poor weather 
conditions the average rms noise in the spectra is above 
0.32 K. With the given sensitivity, for the region of interest 
we do not detect any of the higher excitation lines at 1 mm. 



http://www.iram.fr/IRAMFR/GILDAS 
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Fig. 1. An IRDC filament along the Galactic plane, spanning across from the lower right (18. 45°, -0.20°) to the upper 
left (19.30°, -0.02°). The top panel shows a 8 /mi SPITZER image and the color stretch is chosen to bring the IRDC 
to prominence. The lower panel presents the intensity weighted peak velocity (moment 1) map of the GRS 13 CO data. 
Regions with velocities outside the 43km/s to 51km/s are clipped. The green contours are from ATLASGAL with 
contour levels at 0.15 J y or 3 a, 0.3 Jy, 0.4 Jy , 0.5 Jy, 0.7 Jy, 0.9 Jy, 1.3 Jy, 1.8 Jy, and 2.5 Jy. The yellow ellipses indicate 
the bubbles as given in Simpson et al.| ( |2012 ). 



2.3. The large scale cloud complex in Hi, VGPS 



2.4. A large-scale view of the molecular cloud, GRS 



The ISM is composed primarily of atomic hydrogen. Hi 
(self) absorption features are usua lly spatially correlate d 
to molecular clouds and IRDCs ( Li fc Goldsmith| 2003). 
We used the 21 cm Hi observations of the VLA Galactic 
Plane Survey (VGPS, |Stil et al.||20Q6| > to study the large 
scale distribution of atomic hydrogen. With the VLA in D 
configuration, the best resolution at 21 cm is 45". The inter- 
ferometer data were supplemented with short spacings from 
the NRAO Greenbank telescope (GBT). To improve the 
survey's sensitivity to extended, low surface bright emis- 
sion, the data were reduced to a synthesized beam of l'at 
a velocity resolution of 1.56 km/s. The data cubes are sam- 
pled to a pixel scale of 18"at 0.8 km/s. The calibration was 



made to be consistent with the NVSS survey ( |Condon et al. 
1998). 



To study the molecular cloud component to its full extent 
we selected 13 CO, a (well known) low density gas tracer 
with a critical density ~1.9x 10 3 cm -3 . While it traces the 
same density range as 12 CO, the isotopologue 13 CO has 
lower abundances and is therefore optically thinner. The 
Boston University-F CRAO Galactic Ring Survey (GRS, 



| Jackson et al. [2006] ) provides spectroscopic 13 CO (1-0) ob- 
the~~ 



servations of the Galactic plane for Galactic longitudes 18° 
< 1 < 58°. The observations were performed with the single 
sideband focal plane receiver SEQUOIA, mounted on the 
Five College Radio Astronomy Observatory (FCRAO) 14 m 
telescope in OTF mode. Its angular resolution is 46", and 
the publicly available maps are sampled to a grid of 22". 
The velocity resolution is 0.2 km/s, the pointing accuracy 
is better than 5". 



3 



J. Tackenberg et al., 2012: Triggered/sequential star formation? 



2.5. The stellar component of the complex, UKIDSS, 
GLIMPSE and MIPSGAL 

To examine the stellar content and stellar proper- 
ties, we employed near- and mid-IR surveys. Both the 
Galactic Legac y Infrared Mid-Plan e Survey Extraordinaire 
(GLIMPSE, |Benjamin eriuT1|2003t and the MIPS Galactic 
plane survey (MIPSGAL |Carey et al.|[2QQ9t are SPITZER 
Legacy Programs. The GLIMPSE survey provides maps of 
the inner Galactic plane for -60° < 1 < 60° at 3.8 /im, 
4.5 /im, 5.8 /im, and 8.0 /im. The resolution is 1.7", 1.7", 
1.9", 2.0", respectively. From MIPSGAL we used the 24 /im 
images with a resolution of 6". Additional near-IR JHK 
data have been taken from the UKIDSSDR7PLUS data re- 



lease ( |Lawrence et al. 2007). 



2.6. The peak of the SED, HiGal 

Covering the peak emission of the spectral energy dis- 
tributio n (SEP) of cold dust the HIGAL/HERSCHEL 
survey (Molinari et al. 2010) provides mid-IR observa- 
tions at 70 /im and 160 /im, and at submm wavelength 
250 /im, 350 /im, and 500 /im. Therefore, they connect 
the ATLASGAL submm observations of the cold dust at 
870 /im to the near- and mid-IR observations of the hot 
dust and stellar radiation. From a black body fit to the 
fluxes, we will estimate the temperature of the dust. The 
observations were carried out in parallel mode with a scan- 
ning speed of 60"/s. Since the official, fully reduced HiGal 
fields are not yet availabl e, we red uced the HiGal raw data 
using HIPE (level to l,|Ott|2010D and SCANAMORPHOS 
(level 1 to 2, vl6.0, Roussel |2012p . The angular resolutions 
at the given wavele ngHisare 10.2", 13.6", 2 3.4", 30.3", and 
42.5", respectively flTraficante et al.[|2011). 



2.7. The ionized gas, MAG PIS and SHS 

In addition to thermal dust and line emission, free-free 
emission can be observed, which traces ionized gas. A cer- 
tain amount of ionizing photons per second is required to 
maintain a given amount of gas in ionized state. This allows 
a characterization of the ionizing source. However, espe- 
cially at cm wavelength another significant emission contri- 
bution is coming from the cosmic ray electrons' synchrotron 
radiation. Here, the spectral index can help to distinguish 
which contribution dominates. 

The Multi-Array Galactic Plane Imaging Survey 
(MAGPIS, |Helfand et aT^[2006| mapped the 1.4 GHz, or 
20 cm, continuum of the first Galactic quadrant partially. 
Using the VLA in B, C, and D configurations allows a reso- 
lution of ~ 6", while the extended structures are pres erved 



by including Effelsberg 100 m data ( Reich et al.||1990[ ). The 
ready reduced maps have a pixel scale of 2"/pix. 

While being very susceptible to extinction, the opti- 
cal H a line traces the down cascade of an electron af- 
ter recombining from ioniz ed to atomic hydro gen. The 



SuperCOSMOS H a survey (Parker et al., 2005) digitized 



AAO/UKST narrow-band H a observations of the southern 
Galactic Plane with a resolution of ~ 1" . 



3. Results 

The employment of the various surveys allows us to draw 
a comprehensive large scale picture of the cloud complex 



associated with the IRDC G18.93. In the following we will 
describe the environment, going from large to small scales. 
Finally we will use velocity information to study the spatial 
relations. 

3.1. A massive IRDG in the vicinity of a bubble 

As mentioned in the introduction, a prominent feature 
along the > 50 kpc long filament is the bubble like structure 
indicated by the yellow ellipse in Fig, [l] A second bubble , 
N24 from |Churchwell et al.| ( |2006| ) , |Deharveng et aL] ( |2010 ) 
marked by the lower yellow ellipse in Fig. [lj is close in pro- 
jection as well. However, while the filament has a velocity 
(v^ sr ) around ~45km/s, the second bubble has a distinct 
velocity at around 60km/s and therefore is at a different 
distance. 

A zoom in on the smaller bubble closer to the filament is 
shown in Fig. |2j It opens towards higher Galactic latitudes 
and only little 24 /im emission is visible within the bubble. 
The o pening of the bubble re sembles a typical 'champagne 



flow' (Stahler & Palla, 2005) with lower density material 
at higher latitudes. On the lower inner edge of the rim a 
bright IRAS source is located together with two bright NIR 
sources. 

Directly at the edge of the bubble at lower latitudes, 
several studies identified infrared dark clouds (IRD Cs) and 
massive dense clumps, e.g. SDC G18. 928 -0. 031 (Peretto 
|fc Fuller! [2009] ), [SMC2009] G1 8.93-0.01 jSchuller et al.[ 
2009T^r BGPS G18.926-0.019 ( [Rosolowsky et al.| |2010p . 
ATlhown in Fig. [2| the ATLASGAL survey at 870 /im 
revealed several dense structures along the bubble. The 
largest structure in F ig. [2] is partly foll o wing the IRDC 



filament. In addition, Cyganowski et al. ([2008) found an 



extended green object (EGO) within the rim of the bub- 
ble, EGO G19.01-0.03 which drives a bipolar outflow and 
is characterize d as a genuine fast accreting mas sive young 
stellar object ( Cy ganowski et alT| 2008|[2011b|a[ ). As shown 
in Fig. [3j the main emission peak of the IRDC shows up in 
absorption up to even 160 /im. While absorption of IRDCs 
usually refers to extinction against the Galactic background 
from PAHs and very small grains, it is unclear whether the 
absorption feature at the longer wavelength of 160 /im im- 
plies a lack of emission or simply no grains hot enough 
to emit in that line of sight. However, the part around 
EGO G019.01-0.03 and all connected continuum emission 
towards lower galactic latitudes does not even show 24 /im 
extinction signatures. 

Within and along the bubble, there are few additional 
complexes dense enough to contain dense clumps and be- 
come visible at ATLASGAL sensitivity. However, since we 
want to study the influence of an expanding Hll region upon 
an already existing filament, we will not discuss the contin- 
uum emission in the region that is clearly above the fila- 
ment. Due to its characteristic absorption we will call the IR 
dark structure defined by the ATLASGAL emission G18.93 
(see gray ellipse on Fig. For reasons we will discuss in 
Sec. |3.2| we keep the ATLASGAL emission aligned with 
the galactic latitude axis connected to EGO G019.01-0.03 
separate. 

It is interesting to note that both the massive young 
stellar object EGO G019.01-0.03 and the IRDC filament 
are at the interface of the bubbles shown in Fig. [I] 
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Fig. 2. Three color image of the bubble connected to the long filament, MWP1G018980+00304, composed from the 
SPITZER 3.6 /im, 8 /im, and 24 fim bands in blue, green, and red, respectively. White contours are from ATLASGAL 
(at 0.3 Jy, 0.4 Jy, 0.5 Jy, 0.7 Jy, 0.9 Jy, 1.3 Jy, 1.8 Jy, and 2.5 Jy) for which the numbers denote the CLUMPFIND peaks, 
marked by the red ci rcles. The yellow solid ellipse marks the bubble with interpolated bubble dimensions adopted from 



Simpson et aL|d2012| ). The two arrows mark IRAS 18227-1227, and the massive proto-stellar object EGO G019.01-0.03. 
While the dashed orange boxes indicate the areas that have been mapped in H 13 CO + with the IRAM 30 m telescope, 
the gray dashed ellipse is drawn around the ATLASGAL contours connected to G18.93. 



3.2. Disentangling the spatial relations of G 18.93 and EGO 
G 19.01 

From 13 CO, partly shown in Fig.[l] it cannot be concluded 
whether the dense dust for the region of interest seen in 
ATLASGAL is part of the 45km/s filament component or 
the 60 km/s component connected to the second bubble. To 
address that question and to disentangle the velocity struc- 
ture of the dense molecular gas, we mapped G18.93 and 
the dust emission around EGO G19.01-0.03 in H 13 CO+, a 
dense gas tracer. 

As shown in Fig. [4j the dense gas has two distinct 
velo city components. Consistent with earlier s ingle spec- 



tra (Cyganowski et al.l 12011b, Wienen et al., 2012), the 



IRDC G 18.93 shows emission at velocities between 43 km/s 
and 47 km/s, while the peak of EGO G19.01-0.03 and the 
ATLASGAL emission directly connected to it have veloc- 
ities between 59 km/s and 63 km/s. For both components, 
the dense molecular gas tracer resembles the structure 



defined by the ATLASGAL emission contours very well. 
Therefore, one can attribute the same velocities unambigu- 
ously to the ATLASGAL emission. That is not obvious from 
the lower density 13 CO alone. In addition, it is obvious that 
G18.93 is part of the long filament and the primary bub- 
ble which is connected to it, while EGO G19.01-0.03, a site 
of massive star formation and its connected material below 
this region, is at velocit i es co nsistent with the larger bubble 
fromlChurchwell et al.l (120061). 



In order to constrain the distance to the filament and the 
connecte d bubble, we emp loy the Galactic rotation curve 
given in |Reid et al.| ([£009). Assuming a common velocity 
of 45 km/s, the kinem atic near distance become s 3.6 kpc. 
From both extinction (Ka inulainen et aT| |2Q1 1 ) a nd c on- 
siderations of the Hi absorption feature (see Sec. 3.5 for 
details) we can conclude that the near d istance is appropri- 
ate. Note that Kainulainen et al. ( 2011[ ) give a kinematic 
distance of 4.1 kpc and even larger distances from velocity 
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Fig. 3. A multi- wavelength zoom on G18.93/m, the continuum peak connected to the prominent IR extinction feature. 
The panels are (from left to right in reading order): 3.6 /im, 8.0 /im, 24 /mi, 70 /im, 160 /im, 250 /im, from GLIMPSE, 
MIPSGAL, and HIGAL. The blue circles in the bottom right of each panel indicate the beam sizes, the green contours 
are from ATLASGAL (at 0.3 Jy, 0.4 Jy, 0.5 Jy, 0.7 Jy, 0.9 Jy, 1.3 Jy, 1.8 Jy, and 2.5 Jy). The colors of the so urce marked 
by a red arrow in the first panel a re not consistent with typical colors of young stellar objects (class 0/1) ( jGutermuth 



et al., 2008, Robitaille et al. 2008). 
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Fig. 4. Color coded map of the H 13 CO + emission in [K km/s] for the velocity regime 42km/s to 49km/s on the left, 
and 56 km/s to 65 km/s on the right. The black contours are from ATLASGAL with levels at 0.3 Jy, 0.4 Jy, 0.5 Jy, 0.7 Jy, 
0.9 Jy, 1.3 Jy, 1.8 Jy, and 2.5 Jy. 
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independent methods. Therefore, error propagation result- 
ing in a distance uncertainty of 0.05 kpc underestimates the 
uncertainty. We expect the uncertainty rather to be on the 
order of 0.5 kpc. 

3.3. Details and aspects of the G18.93 complex 

So far we have shown that there is a large filament along 
the Galactic plane that has a consistent velocity structure 
over more than 50 pc. At 1 = 18.93°, b = 0.03°the filament 
becomes very opaque, shows high mid-IR extinction and 
large gas column densities. Part of it shows IR absorption 
up to 160 /im. Just above G18.93, a bubble is visible. Its 
PAH rim visible at 8 /im has a direct interface with G18.93. 
We will show later that the bubble is an expanding Hll 
region. Therefore, G18.93 is an ideal source in which to 
study the influence of massive stars on starless clumps and 
search for imprints of triggering. 

3.3.1. The dust component of G18.93 



As mentioned in Sec. ^2 and shown in Fig. [4j the densest 
gas (and most of the molecular hydrogen) in the complex is 
concentr ated in the regions seen by ATLASGAL at 870 /im 
(see Sec. O). 

In order to identify column density peak positions 
and associate masses to individual clumps, we used the 



CLUMPFIND algorithm of [Williams et 
at a 6 (j contour of 0.3 Jy, we chose 



lumps, we 
FaD([l994 
A0.5TO7 



(1994). Starting 
5T 



0.7, 0.9, 1.3 

1.8, and 2.5 Jy as extraction thresholds. ( For a detailed dis- 
cussio n of the extraction thresholds see [Tackenbe rg et al., 
2012[ ) The positions of the extracted clumps are listed in 
Table [I] and plotted in Fig. [2] We denote the prominent 
absorption dip and most massive clump G18.93/m. 

Assuming optically thin thermal emission, the H2 gas 
column density can be calculated from the dust continuum 



emission via 



N = **± 



(1) 



with the gas-to-dust ratio R= 100, Fa the flux at the given 
wavelength, B^(A,T) the blackbody radiation as a function 
of wavelength and temperature, m# 2 the mass of a hydro- 
gen molecule, and the beam size Q. Assuming typical beam 
averaged volume densities in dense clumps of 10 5 cm -3 and 
dust grains with th in ice mantles, the dust mas s absorp- 
tion coefficient from Ossenkopf & Henning ( 1994 ) becomes 
n = 1.42 cm 2 g _1 at 870 /im. The temperatures are taken 
from the HiGal temperature map described in Sec. 3.3.2 



Because of the sparse resolution, we adopted the temper- 
ature at each peak position. The peak column densities at 
the spatial resolution of ATLASGAL (19.2") are then be- 
tween 1.0 x 10 22 cm -2 to 4.2 x 10 22 cm -2 . 

With the distance as additional parameter and the in- 
tegrated flux of the clump, the mass can be calculated in a 
similar way as given above, 



Rd 2 F; 



X.tot 



gas 



B x (X,T)k 



(2) 



Assuming the kinematic distance of 3.6 kpc, the total mass 
becomes ~ 870 M , with clump masses between ~ 30 M 
and 250 M . Individual column densities and masses are 
listed in Table [l] The given radius is the effective radius 




19.10° 19.00° 18.90° 

Galactic Longitude 




18.98° 18.96° 18.94° 18.92° 

Galactic Longitude 



Fig. 5. Color coded temperature maps in [K] of G18.93, 
shown by the white contours from ATLASGAL (at 0.3 Jy, 
0.4Jy, 0.5Jy, 0.7Jy, 0.9Jy, 1.3 Jy, 1.8 Jy, and 2.5 Jy). The 
beam size is indicated in the top right corner. While the 
top panel presents the large-scale environment, the bottom 
panel zooms into our central IRDC region. 



for equating the pixel area of each clump with a theoretical 
circular area, as calculated by CLUMPFIND. 



3.3.2. Temperature map of G18.93 

The dust temperature is not only a necessary quantity 
for calculating column density and mass estimates, but 
it is an important physical parameter. As described in 
Sec. |2.6| we can use the HiGal/HERSCHEL survey to- 



gether with ATLASGAL data to estimate the dust temper- 
ature. Smoothing all data to the same resolution, we used 
frequ ency-dependent, optically th in dust emission models 



from |Ossenkopf fc Henning 
single temperature Planck 



( 1994 ) to fit every pixel with a 



unction. (See also |Ragan et al] 

L 



|2012[ ) A common problem in the context of HERSCHE 
PACS/SPIRE data is the unknown background contribu 
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Table 1. CLUMPFIND decomposition of dust clumps extracted on ATLASGAL 870 /im map. 



Sub- clump 


Gal. Ion. 


Gal. lat. 


RA(2000) 


DEC(2000) 


HiGal 


Peak column 


Mass 


Angular 


Radius 


name 










peak temp 


density 




radius 






[°] 


[°] 


[ hh:mm:ss.s ] 


[ dd:mm:ss ] 


[K] 


[ 10 22 cm" 2 ] 


[ M ] 


[ "] 


[pc] 


1 


18.9259 


-0.0158 


18:25:32.4 


-12:26:47 


23.0 


4.2 


248 


27 


0.46 


main 


18.9325 


-0.0292 


18:25:36.0 


-12:26:48 


21.7 


3.5 


276 


29 


0.50 


2 


18.9125 


-0.0158 


18:25:30.9 


-12:27:29 


22.3 


2.4 


104 


21 


0.36 


3 


18.9459 


-0.0325 


18:25:38.3 


-12:26:11 


22.5 


2.2 


107 


22 


0.38 


4 


18.9425 


-0.0308 


18:25:37.6 


-12:26:19 


22.4 


2.1 


57 


15 


0.27 


5 


18.9658 


-0.0342 


18:25:41.0 


-12:25:10 


23.0 


1.5 


43 


15 


0.26 



tion. Se veral efforts have b een made to solve tha t prob- 
lem, e.g. |Stutz et"aLl ( |2010[ ) jBattersby etaL] ( |2011[ ). In the 
context of determining the temperature we tested differ- 
ent background levels based on Gaussian fits to the noise 
distribution in regions with no or little signal. While the 
absolute temperatures differ up to 10%, the relative tem- 
perature distribution is very similar. Therefore we refrain 
from subtracting any background and focus on the relative 
changes. 

Furthermore, we have compared temperature maps us- 
ing both HiGal and ATLASGAL, smoothed to a common 
resolution of 37", to temperature maps with only three 
bands, 70 /im, 160 /im, and 250 pm at a resolution of 19". 
Although it is ambitious to fit a curve to only three data 
points, the peak of the SED is covered and can be recon- 
structed. The ATLASGAL data has been omitted since 
the missing background subtraction implies problems with 
the calibration and worsens the reconstruction of the peak. 
Comparing the temperatures of both SED fits on regions 
where differences due to the beam size are negligible, we 
find very good agreements. We want to point out that omit- 
ting the longer wavelength data biases the absolute results 
towards higher temperatures. However, the relative distri- 
bution in the temperature maps are preserved. Therefore, 
Fig. [5] shows the temperature maps at the better resolu- 
tion of 19". To have more reliable absolute temperatures 
for deriving column densities and masses, we use the tem- 
peratures derived on all HiGal plus the ATLASGAL data. 

As displayed in the top panel of Fig. [5j the dust around 
the IRAS source is hottest and produces a large scale tem- 
perature gradient. Nevertheless, all ATLASGAL peaks are 
colder than their neighborhood, and the temperature dip 
towards G18.93/m is the largest. 



3.3.3. G18.93/m: a starless clump 

Despite its p rominence as absorption feature at SPITZER 
wavelengths, Tackenberg et al. (2012) did not list G18.93/m 
as a starless clump. Because of a peak in the extended 24 /im 
emission within the clump boundaries, the clump had been 
rejected for consistency. In this detailed study however, we 
refrain from such rigorous methods and allow an individual 
inspection of the clump. Using the GLIMPSE color crite- 



quantify our sensitivity we estimated the luminosity of ob- 
jects that still could be embedded in the d ust. Using the 
same assumptions as explained in Sec. |3.3.2[ we here fitted 
two blackbody functions to the point source detection lim- 
i ts of the 24 am an d 70 /im images of 2mJy for MIPSGAL 
( |Carey et al.] [20091 ) and ~ 10 mJ y for the 70 A™ HiGal field 
together with the flux estimates for the absorption peak 
of G18.93/m at 160 /im and 250 /im, 52.0 Jy and 45.1 Jy. 
The warmer component then has only 0.1 L , correspond- 
ing to a low-mass star of 0.14 M . This can be considered as 
an upper limit for any potential embedded and undetected 
source. 

Additionally indicative of true starless clumps is the ab- 
sence of SiO emission. Since SiO traces outflows it is com- 
monly used to differe ntiate between starl e ss clumps and 
star forming clumps (iMotte et al. 2007, Russeil et al 



2010). For the entire G18.93 complex we find no SiO down 
to column densities of ~ 8 x 10 11 cm -2 . 

In order to understand whether G18.93/m is pre-stellar 
or only a transient object, we compare its virial mass to 
its dust mass. While pre-stellar clumps are gravitationally 
bound and therefore will eventually form stars in the future, 
transient clumps are not gravitationally bound. Within 
such objects, smaller (here unresolved) fragments can still 
collapse, but the c lumps itself will drift apart. As we will 
discuss in Sec. 4.5 , we find a double peaked H 13 CO + profile 
at the ATLASGAL peak position of G18.93/m. As shown in 
Fig. [TT] the line width of the component we attri bute to the 
IRDC is £v = 2.1km/s. (For details see Sec. [45l) To calcu- 
late the virial mass we use the equation given in MacLaren 
et al.| ( |1988| ) , M vir = k R 5v 2 . For the source radius R we use 
0.5 pc as given in Table [I] and the geometrical parameter k 
depends on the density structure of the clump and is k = 190 
for p oc 1/r, k = 126 for p ex 1/r 2 . The mass then becomes 
M vir = 420 M for a geometrical parameter k = 19 or p oc 
1/r, a nd M vir = 28 with k= 126 or p oc 1/ r 2 . |Beuther 



et al.|fl2002D and |Hatchell fc van der T^k| (|2003) find typical 
density distributions in sites of massive star formation of 
p oc i a with a~-1.6, in between both parameters. Doing 
a linear interpolation, the suggested virial mass becomes 
360 M . Considering all uncertainties, a comparison to the 
estimated dust mass for G18.93/m of 280 M does not al- 
low a conclusive distinction whether clump G18.93/m is 
gravitationally bound or not. However, as discussed in Sec. 



ria given in Gutermuth et al. (2008), none of the cataloged 3.3.2 the dust temperatures derived from SED fitting to 



near- and mid-IR sources within the clump is classified as 
young stellar object (class 0/1). Us ing similar GLIMPSE 
color criteria, Robitaille et al. (2008) did not identify young 
stellar objects projected onto the clump either. For even 
younger sources not yet visible at near-IR wavelength, Fig. 
[3] shows that there is no point source longwards of 8 pm 
close to the ATLASGAL peak of G18.93/m. To further 



HERSCHEL data are very uncertain due to missing back- 
ground levels. Co mparing gas temper atures measured with 
NH 3 observations (Pillai et al. , 2006 ) to dust temperat ures 
from HERSCHEL observations (|Henning et al. , 2010), we 



find that despite the efficient dust cooling the ammonia 
temperatures are usually lower. Thus, the given temper- 
atures are only upper limits and the clump masses may 
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19.05° 19.00° 18.95° 18.90° 

Galactic Longitude 

Fig. 6. The ionized gas within and around the bubble. 
On top of the GLIMPSE 8/im image, the red contours 
show the H a emission from the SuperCOSMOS survey, the 
blue contours represent the MAGPIS 20 cm emission, and 
white contours show the ATLASGAL emission (at 0.3 Jy, 
0.4 Jy, 0.5 Jy, 0.7 Jy, 0.9 Jy, 1.3 Jy, 1.8 Jy, and 2.5 Jy). The 
contour levels for the MAGPIS data are 0.002 Jy/beam, 
0.003 Jy/beam, 0.004 Jy/beam, and 0.005 Jy/beam. The ar- 
row indicates the peak position of the Hll region identified 
by |Lockman| ( |1989D . 



be higher. Therefore, the virial analysis is consistent with 
G18.93/m being gravitationally bound, and pre-stellar. 

3.4. The bubble 



According to |Simpson et al. ( |2012| ), the bub ble c on- 
filament 



3.1 



nected to the filament and described in Sec. 
MWP1G018980+000304, has an effective radius of 3.44'. 
However, as mentioned before, its appearance is more sim- 
ilar to a 'champagne flow' (cf. Fig. |2l). Since the expansion 
velocity of an Hll region depends on the surrounding den- 
sity, an inhomogeneous medium leads to asymmetric bub- 
bles. If the ionizing gas reaches the edge of its parent molec- 
ular cloud or enc ounters dense gas, it w ill fan out towards 
the empty space (iStahler fc Palla| [2005 ) . Therefore, the ef- 
fective radius and the major- and minor-axis can describe 
the current appearance of a bubble, but cannot describe 
the expansion history. 

Thus we have to ask: What is the origin of 
MWP1G018980+00304? Galactic bubbles are usually pro- 
duced by OB stars or clusters, driving an expanding Hll 
region, or by stell ar winds and radiation pressure of late 
B-type stars ([E lmegre en fc Lada[ 19771 Churchwell et al.[ 
20061 Peharveng et al.| |2010[ jSimpson et al.| |2012p 7Tbr 
both processes, the exciting source is not necessarily in the 
middle of the bubble. 



Often, (ultra) compact Hll regions (UCHlls) are found 
connected to the bubb le driving source, but e.g. CORNISH 
2008[ ) does not find an UCHn region in- 
However, despite its high sensitivity at 



(Purcell et al. 
side the bubb 



6 cm of better than 2 mJy, due to interferometry tech- 
niques, _CORNISHis not sensitive to sources larger than 
12 // ( [Hoare et al.[|2012| ). Searching for more extended struc- 
tureTTLockman (1989) conducted a H a radio recombina- 
tion line single dish survey at 3 cm and detected an Hll 
region at 1=18.954, b=-0.019. As shown in Fig. [6j this is 
right between the IRAS source and G18.93, but measured 
with a beam of 3' . Although they measure a source ve- 
locity of vlsr = (52.3±1.7) km/s, we believe that the bub- 
ble and the Hll gas are spatially connected. Indeed, veloc- 
i ty offsets between the ionized gas and CO are common 
([Blitz et all |1982[ |Fich et al.[ |1982[ ). For the same Hll re- 
gion IKuchar & Clarkl 01997 measure the flux at 4.85 GHz 
with a beam width of 4.2' . For an estimated source diam- 
eter of 7.0 / they fin d 2146 mJy. Using the formula given in 
Kurtz et al. (19941) and assuming a correction factor for the 



optical depth of a = 0.9922 (corresponding to a fr equency 
v = 5 GHz, and tem perature T =9000K as given in Mezger 



& Henderson 1967), this flux translates to the number of 
Lyman continuum photons n^ =log(N^) =48.4s _1 . 

The blue contours in Fig.|6]show the MAGPIS GPS data 
at 20 cm. The emission peak of the high resolution data at 
20 cm agr e es wel l with the extrapolated peak position from 



Lockman (1989) in the sense that it also lies between the 
IRAS source and G18.93. However, the GPS data shows 
20 cm emission towards the IRDC and beyond. The higher 
latitude flux within the contour level 3 x 10 -3 Jy/beam, 
without the tail towards lo wer latitudes, becomes 1.5 Jy . 
With a = 0.9951 at 9000K ( |Mezger fc Henderson) |l967| , 
that converts to a Lyman continuum flux of n^ = 48.2 s . 
A slightly larger contour drawn along the visual extent of 
the emission provides a flux of 2.5Jy, thus 2/3rd bigger. 
However, with n^ = 48.4 s _1 the logarithm does not change 
as much and the given differences can be considered as un- 
certainty in our mea surements. Therefore, the results of 
|Kuchar fc Clark| ( |1997[ ) and our measurements agree within 
the uncertainties. 

As mentioned in Sec. |2.7[ the synchrotron radiation 
of cosmic ray electrons can also contribute to the cm 
continuum emission. However, the spectral index a (with 
§ v (x.v~ a ) of the thermal free- free emission is negative in 
this regime, while the spectral index of the synchrotron ra- 
diation is positive. Therefore, we can use the two indepen- 
dent measurements of the continuum emission to at least 
determine the sign of the spectral index. The slope between 
the 20 cm, or 1.5 GHz, data point and the 6 cm, or 4.85 GHz, 
data point is positive, suggesting a negative spectral in- 
dex, which is in agreement with the expected spectral in- 
dex for thermal free-free emission. Therefore, the emission 
is dominated by free-free radiation and we can neglect any 
synchrotron contribution. Nevertheless, with only two data 
points we cannot distinguish whether the assumption of op- 
tically thin emission is correct. Therefore, the given Lyman 
continuum fluxes are lower limits. 

While the cm continuum is more prominent towards the 
dust continuum emission, large parts of the bubble are filled 
by SuperCOSMOS H a emission, see Fig. [6j As mentioned 
in Sec. |2.7[ H a is a complementary tracer of ionized gas. 
The difference in the spatial distribution can be explained 
by the optical depth. Since the optical R a line is suscep- 
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tible to extinction, R a can only be detected in regions of 
low visual extinction. However, the near-IR extinction map 
shows that the visual extinction towards the ATLAS GAL 
emission peaks are larger than Ay = 25 mag. 

Using the averaged number of Lyman photons, 
iiLy = 48.3 s -1 to determine the spectral type, the exciting 
source needs to be at least a main seq uence star of type 08.5 
( |Panagia] [l973, M artins et al.| [2005). The extended appear- 
ance (and the connected non-detection within CORNISH) 
and the morphology of the Hll suggests that the ionizing 
source is not extremely young. Nevertheless, the diameter of 
the bubble does not allow for extremely old bubble-driving 
sources. 

From the mid infrared data, a good candidate for the 
ionizing source is IRAS 18227-1227, marked in Fig. |2j It is 
very bright at 24 /im, but by 3.6 jam it is no longer the 
brightest source within its neighborhood. At higher lat- 
itudes, - 16" above IRAS 18227-1227, a near IR source, 
hardly visible at 24 /im, is another potential candidate to 
drive the bubble. Although the additional source appears 
blue in the four IRAC bands from GLIM PSE, the color 
criteria given in |Gutermuth et ah] ( |2008| ) identify it as 
young stellar object. At the even shorter wavelengths of 
the 2MASS and UKIDSS near IR JHK surveys, the second 
source has a bright neighbor. Nevertheless, its missing de- 
tection in the GLIMPSE survey and its blue appearance at 
NIR wavelength indicates that it is evolved. However, the 
near infrared colors do not allow a classification or mass de- 
termination. While for one its classification as young stel- 
lar object (class I) hints to non photospheric emission at 
near infrared wavelength, the second source shows a sig- 
nificant color excess in the JHK color-color diagram indi- 
cating non photospheric emission as well. Therefore we can 
not identify the exact ionizing source with the current data. 
Nevertheless, as shown before from cm free-free emission, 
it has to be of spectral type 08.5 or earlier. 

Another result one can determine from the different 
tracers of the ionized gas are the densities above and below 
the ionizing source. While we have material that is opti- 
cally thick to visible light suggesting dense gas towards the 
dust continuum emission, we have optically thin low density 
material within the bubble and above it. That explains the 
position of the ionizing source and the 'champagne flow'. 

3.5. The photon dominated region: a layered structure 

What is the observable impact of the Hll region on the 
dense filament, described in Sec. [3]? 

During the evolution of an Hll region the temperature 
difference between the hot ionized gas and the cold environ- 
ment drives a supersonic shock front beyond the ionization 
front. In addition, at the interface between an Hll region 
and the neighboring gas a photon dominated region builds 
up. While hydrogen ionizing radiation with energies above 
13.6 eV produces the Hll, beyond the ionization front other 
molecules with lower ionization energies may still be ion- 
ized. However, connected to the photon domin ated region 
one expects a lay ered structure of Hll - Hi - H2 ( Hollenbachj 
fc Tielens 1997). For our bubble we can directly observe 
the layers between the presumably ionizing source and the 
dense gas of G18.93/m. 

While the ionized gas and its structure is directly traced 
by the cm free- free emission, the ATLASGAL emission con- 
nected to G18.93/m represents the cold dust, which is em- 




Velocity (k 



Fig. 7. HI VGPS spectrum of G18.93/m in black with the 
GRS 13 CO spectrum in red on top. The strong absorption 
feature of G18.93/m at ~ 45 km/s in the Hi spectrum agrees 
very well with the 13 CO signature. The smooth line indi- 
cates the 3rd order fit to the Hi spectrum that has been 
used to subtract the 'continuum'. (For details see text.) 



bedded in dense H2. In order to trace the morphology of 
the atomic hydrogen we employ the VGPS Hi survey. 

(Hot) Hi is so abundant in the Galaxy, that along the 
Galactic plane, background Hi emission is present at all ve- 
locities. Therefore, cold sources in the foreground appear 
as absorption features superimposed on the background Hi 
emission. However, since the Hi distribution and temper- 
ature vary strongly within the Galactic disk, the observed 
background varies strongly as well. In addition, individual 
Hi clouds often have sufficiently high column densities to 
become self absorbing. This complicates the identification 
of complexes in the Hi data. 

However, as shown in Fig. [7J the absorption dip for the 
dense gas around G18.93 is very prominent. In order to 
study the distribution of the atomic hydrogen, we quanti- 
fied the absorption feature by the following process: (1) We 
fit the background in the spectra with a third order polyno- 
mial omitting the velocity range from 40 km/s to 50 km/s. 
This determines the large scale background variations, but 
preserves the signals on small scales. (2) We subtract the 
fitted function as a baseline. (3) We invert the spectra and 
fit the now emission peak with a Gaussian. The integrated 
area of the Gaussian for a constant temperature is pro- 
portional to the column density. (NOTE: We stress that 
the authors do not want to convert the so quantified mea- 
sure of Hi into physical meaningful units, but only take it 
as relative value.) Doing this procedure on all Hi spectra 
around G18.93/m, we obtain the qualitative distribution of 
the atomic hydrogen. 

A problem with this method arises from the free-free con- 
tinuum. Due to the additional background emission towards 
the Hll region, the absorption dip becomes stronger for the 
same Hi column densities. Therefore, the interpretation of 
the Hi distribution needs to be considered with caution. 
However, since the peak of the Hi and Hll emission do not 
agree, the suggested Hi distribution seems to be dominated 
by the Hi column density. 

Figure [8] shows three physical states of hydrogen. 
Between the ionizing source and G18.93 a layered struc- 
ture is visible, with first the ionized gas, then the atomic 
hydrogen, and then the molecular gas within the IRDC. 
However, one should keep in mind that the beam of the 
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Fig. 8. A layered structure of the hydrogen phase. On top 
of the 8 /im GLIMPSE image in gray, the blue contours rep- 
resent the MAGPIS GPS 20 cm data, the red contours show 
a measure of the atomic hydrogen at velocities of G18.93 
from the VGPS 21cm line (for details see the text), and 
the green contours are the cold dust from the ATLAS GAL 
survey, hence the molecular hydrogen. The dashed lines are 
to guide the eye and indicate the peak positions of the dif- 
ferent data. The orange arrow points to 8 /im excess emis- 
sion at the edge of the dense gas. The top left circle in- 
dicates the beam size of the VGPS data, while the bot- 
tom circle shows the beam of ATLASGAL. In this figure, 
MAGPIS is smoothed to the resolution of ATLASGAL. The 
contour levels for the GPS 20 cm data are 2mJy/beam, 
3mJy/beam, 4mJy/beam, and 5mJy/beam, and for the 
ATLASGAL 870 /im emission 0.3 Jy, 0.4 Jy, 0.5 Jy, 0.7 Jy, 
0.9 Jy, 1.3 Jy, 1.8 Jy, and 2.5 Jy. 



Hi data is on the same order as the spacing between the 
different peaks. 



In the following Section we will explore whether we can 
see any imprints of the Hll region on the starless clump 
G18.93/m. 



4.1. Comparison of G18.93/m with typical high-mass starless 
clumps 



As discussed in Sec. 3.3.3 , G18.93/m is a starless clump 
with no IR emission up to 160 jam. While the virial anal- 
ysis does not allow a firm conclusion as to whether it is 
bound or not, but if it is bound it would be a proto-type 
pre- stellar clump. In terms of mass and size it does not 
st and out. If we compare i t to the starless clumps found 
in |Tackenberg et aL (2012), its clump mass of 280 M is 
comparable to the average mass of 315 M they find. Note 
th at with the lower dust o pacity of ft = 0.77 cm 2 g _1 used 
in Tackenberg et al. ( |2012[ ) as well as the low temperature 
of 15K they assume for starless clumps, the clump mass 
becomes 900 M . Nevertheless, for peculiar sources such as 
G18.93/m, the elevated dust temperature needs to be taken 
into account. 

If we assume a spherical clump we can calculate the av- 
erage volume density of G18.93/m to be 1 x 10 4 cm -3 . This 
is significantly smaller than the 5. x 10 4 cm" 3 for the sam- 



ple of starless clumps presented in Tackenberg et al. (2012). 
However, from the detection of H 13 CO^ with a critical den- 
sity of ~ 1.8 x 10 5 cm -3 , we know that at least in the cen- 
tral regions the density needs to be higher. Nevertheless, 
since the average densities of the starless clumps have 
been calculated with the same assumptions, the results 
should be comparable. Fo r the formation of high-mass 
stars, Krumholz & McKee (2008) require a peak column 



density of 3 x 10 z,3 cm _z . With a peak column density of 
4 x 10 22 cm -2 , G18.93/m does not meet their requirements. 
However, the measured column densities are beam averaged 
values, true peak column densities are expected to be higher 
( |Vasyunina et aL 2009). Although the expanding Hll re- 
gion might have influenced the starless clump G18.93/m, it 
does not differ from other starless clumps. 



4. Imprints of triggering? Discussing the 
interaction between G18.93/m and the 
expanding Hn region 

So far, most studies of triggered star f ormation (e.g. 



Deharveng et al.||2003 Zavagno et al.| |2006) have searched 
for young stellar objects connected to bubbles and Hi I re- 
gions. To understand the statistical significance of trig- 
gered star formation, current and future studies make use 
of the rising number of Galactic plane surveys i Z avagno 
et al.[|2010[ Deharveng et al.||2010| |Kendrew et al.[|2012p . 
Although it has been shown in the past that not all stars on 
the border of Hll regions are necessarily formed by trigger- 
ing, it becomes clear that triggering might have a significant 
effect on star formation. 

As a next step towards a better understanding of trig- 
gered star formation we need to observationally identify the 
mechanisms that govern the interaction between expanding 
Hll regions and the molecular gas. In this context, detailed 
studies of individual regions are needed to rev eal signatures 
of tri ggering before star formation sets in (e.g. Bieging et al. 
|2009| . 



4.2. The photon dominated region: a layered structure 

While the general distribution of the atomic hydrogen fol- 
lows the ATLASGAL emission, the peak of Hi is elongated 
towards the gap between the ionizing source and the molec- 
ular hydrogen, shown by Fig. [8] Towards G18.93/5 (cf. Fig. 
|2j, the 20 cm free-free emission is no longer parallel to the 
cold dust, but crosses the IRDC. However, both extinction 
and dust emission are significantly lower at the intersec- 
tion. That implies that there is less dense gas which allows 
the ionized gas to more easily escape beyond G 18. 93/5; the 
ionizing radiation is simply no longer blocked across the full 
height of the filament. 



In the context of star formation, |Glover fc Clark (2012) 
have shown that the composition of the gas has hardly 
any influence on the star formation efficiency. They suggest 
that molecular gas is not a prerequisite for the formation 
of cold clumps/cores, but that the high (column) densities 
required also prohibit the destruction of molecules. Hence 
molecular gas is not a requirement, but form simultane- 
ously. Therefore, the layered structure does not affect the 
star formation within G18.93. Instead, shielding of the ra- 
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diation field is most important and the layered structure 
proves the effectiveness of the shielding. 



4.3. The temperature distribution 



As described in Sec. |3.3.2[ the hot ionizing source as well 
as the IRAS source heat up their environment and produce 
a temperature gradient across G18.93 (cf. Fig.j5|. On large 
scales shown by the top panel of Fig. [5j beyonclthe heating 
from the ionizing source the fitted dust temperatur e drops 
to values be tween 19 K and 20 K, typical for the ISM (|Reach| 



et al. 1995). Within, the temperatures towards the contin- 
uum peaks stand out against their direct vicinity. However, 
all dense clumps are warmer than the general dust/ISM 
temperature. 

A different visualization of the temperature gradient be- 
comes visible if we compare the mid-IR wavelengths. Figure 
[9] shows a cut through G18.93/m for different wavelengths, 
ranging from 8 /im to 870 /im. While the longest wave- 
lengths trace the column density peak of the cold dust, the 
SED peak moves along the temperature gradient towards 
the heating source. Therefore, the shorter wavelengths are 
offset towards the heating source as well and we can trace 
its transition directly if we employ the full resolution at 
all wavelengths. Only the 8 jam band is not consistent with 
that picture. Since it is dominated by PAH emission, its 
peak is influenced by the ionizing radiation. 

IRDC temperatures elevated above the general ISM 
temperature of ~ 18 K are different from t ypical regions 
of both low- and high-mass star formation ( |Peretto et al.| 
[20101 jBattersby et al.|[20TTj Nielbock et al., submitted, 
Launhardt et al., in prep.) in which their temperatures (or 
the ones of dense cores) drop to temperat ures below the 
ISM v alue. In contrast to the above studies, |Beuther et al.| 
d2012| ) find values for IRDC 18454 similar to what we find. 
In IRDC 18454, the mini starburst cluster W43 raises the 
general temperature of the IR DC complex to simila r tem- 
peratures as we find for G18.93. Beuther et al. (2012) point 



out that the elevation of the dust temperature raises the 
Jeans length/mass, and therefore speculate that high- mass 
star formation may be favored. Hence, it could be that it 
is not the expansion of the Hll region and the connected 
shock front that promotes the formation of OB associa- 
tions on the rims of bubbles, but the effects of heating on 
the environment. 

Therefore we conclude that the bubble does have an 
effect on the temperature of the IRDC, but we cannot dis- 
tinguish whether this is only due to the heating source or 
also because of the shock front. 



4.4. Does the H11 region change the shape of the IRDC? 



Elmegreen ( 1998| ) explains the formation of fingers (or ele- 
phant trunk like structures) as a shock wave runs over den- 
sity enhancements and this has been directly observed (e.g. 
|Motte et aLl|2010[ ). For the interface between G18.93 and 
the neighboring bubble no such structures are visible in Fig. 
[2] nor at other wavelengths available. 

In the picture of the 'collect and collapse' model the 
Hll region pushes material along. In the vicinity of a pre- 
existing IRDC, the shock wave mainly penetrates one side 
of the IRDC, the one facing the Hll region. Thus the IRDC 
should become asymmetric, with a steeper density profile 



towards the bubble. Such an asymmetry should be reflected 
in the projected profiles of the IRDC. 

Figure 10 shows profiles of various gas tracers along 
two lines. While the profile through clump G 18. 93/1 shows 
a second peak along the profile and is therefore unsuited 
for comparing the shape of the profile, the cut through 
G18.93/m seems to be well suited. Figure 10 shows that 
13 CO is more of an envelope around the entire complex 
and the Hll peak is enclosed within the low density gas. 
Different from that, the dense gas peak from ATLAS GAL 
data is clearly narrower. However, fitting both sides of the 
sub-mm continuum intensity distribution of G18.93/m with 
an exponential function independently, we find no difference 
between both sides of the profile. At the resolution of the 
IRAM H 13 CO + maps we do not have sufficient data points 
to do statistically meaningful fits to both slopes, but on the 
sparse data available, again we find no difference. 

Therefore we conclude that at least for G18.93/m the 
bubble does not seem to influence the density profile of the 
IRDC. 



4.5. Imprints of the shock in the dense gas? 

After the fast (super sonic) initial expansion of an Hll re- 
gion to the size of a Stromgen sphere, the shock's velocity 
quickly drops below the sound speed vs of the medium. 
Typical values of v^ are close to lOkm/s for the hot ion- 
ized gas. As the shock front hits the dense gas of G18.93, 
the expansion velocity drops dramatically, while the lower 
density gas towards higher latitudes allows an accelerated 
expansion producing the 'champagne flow'. However, even 
the slow shock front should still leave its imprint in the 
dense gas. 

In order to search the potentially shocked gas for veloc- 
ity peculiarities, we examined H 13 CO + spectra at various 
positions. 



For three clumps, Fig. 11 compares H 13 CO + spectra 
within the 8 jam rim to spectra at the clump's peak position 
and positions away from the Hll region. Closest to the shock 
front, the spectra of G18.93/m and G 18. 93/1, the two right 
panels of Fig. [TT] consist clearly of two distinct velocity 
components of similar strength. Since H 13 CO + is expected 
to be optically thin, we can exclude self absorption. The 
spectra at the position of the ATLASGAL continuum peaks 
show two components at the same velocities. However, here 
one component is significantly stronger than before. At the 
third position, away from the Hll region, G18.93/m shows 
only a single component, while we do not detect signal for 
G18.93/1. 

A possible interpretation of those spectra could be that 
one component traces the IRDC, while the second compo- 
nent is shock induced. In that context the rising component 
towards the peak of the IRDC traces the density distribu- 
tion of the dense gas. Its distribution is almost symmet- 
ric perpendicular to the filament. The second component 
is bright towards the Hll region, but due to the IRDC's 
high density, the shock has not propagated beyond the peak 
and therefore no shock component is visible away from the 
shock front. For the spectra taken across G18.93/4, within 
our sensitivity limits we cannot find similar clear imprints. 
Although the spectrum at its peak position has no pure 
Gaussian profile and a second component is required to 
fit it, a similar interpretation as above does not work. In 
Sec. |4.2[ we argued that the layered structure breaks down 
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Fig. 9. Profiles of mid IR data along the cut indicated by the yellow arrow in the left panel. While the center of the 
cut in the right panel is G18.93/m, its upper end is chosen to be towards the potential ionizing sources. The different 
colored profiles in the right panel are of increasing wavelength, starting with 8 fim at the top going to 870 fim at the 
bottom. The intensities are not to scale, but adopted to fit in the panel. The green vertical line indicates the position of 
the ATLASGAL continuum peak, the yellow arrow corresponds in direction and length to the arrow in the left panel, 
in which the background image is a GLIMPSE 8/im image with ATLASGAL contours (at 0.3 Jy, 0.4 Jy, 0.5 Jy, 0.7 Jy, 
0.9 Jy, 1.3 Jy, 1.8 Jy, and 2.5 Jy) on top. 




Galactic longitude / [ ° ] projection to Gal. lat. / [ ° ] 

Fig. 10. Profile of 13 CO (black), ATLASGAL (red), and H 13 CO + (dashed blue) along the two yellow lines shown on top 
of the 870 {im map in the left plot. The position of the peak of Hll emission is indicated within both profiles (right) and 
the map (left). The left panel also gives the two main IRDC peaks. The blue asterisks indicate the position of the IRAS 
source (lower position) and UKIDSS source (higher position). The green contours from ATLASGAL are at 0.3 Jy, 0.4 Jy, 
0.5 Jy, 0.7 Jy, 0.9 Jy, 1.3 Jy, 1.8 Jy, and 2.5 Jy. 



clumps at higher longitudes, because of their lower density. 
For the same reason, the shock might have had less impact 
on the dense gas at the center of the filament. Therefore, 
the imprints in the spectra are not as clear. 



It is interesting to note that Klessen et al. (2005) finds 
similar imprints on quiescent cores. For their modeling of 



dense cores in which the turbulence is driven on large scales, 
the velocity dispersion a tur b of their simulated cores show 
bow-shaped enhancements (rims) around the continuum 
emission. In our picture, instead of a large scale convergent 
flow, the shock introduces tur bulence on large scales. The 
fact that Klessen et al. ( 2005[ ) measure an increase of line- 
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Fig. 11. H 13 CO + spectra at positions indicated by bold asterisks in the top panel. The positions are chosen to lie along 
the connection between the UKIDSS near-IR source and the ATLASGAL peak position, with the first spectra at the 
peak of the GLIMPSE 8/im emission. From left to right the clumps are G18.93/1, G18.93/m, and G18.93/4. The upper 
spectra are at positions closer to the ionizing source. Where possible we fitted one, two, or three gaussian components 
to the spectra. The fitted line widths and peak positions are given in the spectra. 



width instead of two separate velocity components could 
well be a result of their method for determining (Jturb- 

However, the Hll region does influence the dense gas of 
the IRDC. The additional velocity component might give 
the IRDC additional turbulent support to build more mas- 
sive fragments. 

Another possible imprint of a shock on dense gas could 
be a broadened line width. Within the limits of our H 13 CO + 
maps, we do not observe such a broadening. 



5. Conclusion 

In the middle of a > 55 pc long filament is G18.93, a promi- 
nent IRDC. While it has the luminous protostellar object 
EGO G19.01-0.03 at one end, the other end is dark up 
to 160 fim. Particularly interesting about G18.93 is its en- 
vironment; it is located at the projected interface of two 
IR bubbles. However, it is only spatially coincident with 
the bubble above the filament, while the bubble at lower 



latitudes has a different velocity. Using H 13 CO + as dense 
gas tracer we unambiguously attribute kinematic distances 
to the ATLASGAL continuum emission. Although very 
close in projection, the high-mass star forming region EGO 
G19.01-0.03 and its cold dust is not connected to G18.93, 
but exists at a different distance. From a Galactic rotation 
curve we determined IRDC G18.93's distance to be 3.6 kpc. 

We used CLUMPFIND to decompose the dense gas seen 
by ATLASGAL into 6 clumps. Together with temperature 
estimates from HiGal SED fitting, we calculated both the 
column density and masses of these clumps. The most mas- 
sive clump is G18.93/m with a gas mass of ~ 280 M . A 
comparison of the gas mass to its virial mass of 360 M 
shows that within the uncertainties it could be gravitation- 
ally bound. To determine its evolutionary stage, we first 
searched the GLIMPSE catalog for young stellar objects 
using IRAC color criteria, but no young stellar object was 
found within the boundaries of G18.93/m. Next we visually 
inspected the longer wavelength images MIPS GAL 24 /im 
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and HiGal and find no point source at wavelength up to 
160 fim. 

In addition, the absence of SiO emission is strongly in- 
dicative of no star formation activity. Therefore, we identify 
G18.93/m as a potential high-mass pre-stellar clump. 

Expanding the SED fitting from single positions to all 
pixels, we produced temperature maps. The sources embed- 
ded in the dust connected to IRAS18227-1227 heat the dust 
and produce a strong gradient across the IRDC. Therefore, 
all IR dark continuum peaks have temperatures above the 
general ISM value of ~ 19 K. 

In context of the bubble we show that the 'champagne 
flow' structure is produced by an expanding Hll region. 
From measuring the radio continuum flux we estimate the 
number of Lyman continuum photons and constrain the 
ionizing source to be at least an 08.5 star. 

From VGPS Hi spectra we extract a measure of the 
Hi column density. Together with the cm continuum data 
and the ATLAS GAL data we identify a layered structure 
between the ionizing source and the IRDC, from ionized 
through atomic to molecular hydrogen. 

Finally, we discuss the imprints of the expanding Hll 
region on the starless clump G18.93/m. Our main results 
are: 

— layered structure: The ionizing source produces a layered 
structure towards the IRDC and therefore changes the 
composition of the hydrogen phase. 

— temperature distribution: The IRAS source and/or the 
ionizing source produce a strong temperature gradi- 
ent across the IRDC. Although the dust temperature 
at the continuum peaks is reduced compared to their 
surrounding, their absolute temperatures are above the 
ISM values. This is atypical co mpared to most othe r 
IRDCs. However, as discussed in |Beuther et al. ( 2012[ ), 
the additional support raises the Jeans length which 
might favor high- mass star formation. 

— density profiles: In the picture of the 'collapse and col- 
lect' scenario we would expect the shock front to steepen 
the density profile of the dense gas. Therefore we com- 
pare the emission profile of the cold gas towards and 
away from G18.93/m. Fitting both wings, we find no 
difference. 

— imprints within the dense gas: Looking at H 13 CO + spec- 
tra across the filament, we find emission that resembles 
the brightness of the ATLASGAL dust continuum. In 
addition we find a second H 13 CO + component, brightest 
within a PAH rim towards the bubble and not detected 
away from the bubble. We speculate that this might be 
a direct imprint of the shock front onto the IRDC. 

Therefore, while the additional heating and the shock 
induced velocity component favor high-mass star forma- 
tion, we do not find evidence for collapse triggered by the 
expanding Hll region. 
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